Quantum Oscillations of Elastic Modiuli and Softening of Phonon Modes in Metals 
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In this paper we present a theoretical analysis of the effect of magnetostriction on quantum 
oscillations of elastic constants in metals under strong magnetic fields. It is shown that at low 
temperatures a significant softening of some acoustic modes could occur near peaks of quantum 
oscillations of the electron density of states (DOS) at the Fermi surface (FS). This effect is caused 
by the Condon magnetic instability, and it can give rise to a lattice instability. We show that the 
most favorable conditions for this instability to be revealed occur in metals whose Fermi surfaces 
include nearly cylindrical segments. 
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Experimental data concerning quantum oscillations 
in various observables in metals under strong magnetic 
fields were repeatedly used in studies of their electron 
characteristics. At low temperatures when the param- 
eter 9 is small compared to unity (8 — 2TT 2 T/hft; T is 
the temperature expressed in units of energy; HQ is the 
value of a cyclotron quantum) these oscillations can ex- 
hibit a rich structure. The latter could be analyzed even 
within a simple isotropic model. For instance, it is known 
that the longitudinal magnetic susceptibility of a metal 
can reveal divergencies at peaks of quantum oscillations 
at low temperatures [1-3]. These divergencies crucially 
depend on interaction among conduction electrons, and 
they indicate a possibility of a diamagnetic phase transi- 
tion in a metal, producing Condon domain structure near 
the oscillation peaks [4]. Also, it was proposed in some 
earlier works that the above features in the magnetic sus- 
ceptibility could give rise to a sensible reduction of elas- 
tic moduli in metals [1,4-6]. Here, we demonstrate that 
the effect of softening of phonon modes near the peaks 
of magnetic quantum oscillations of the electron DOS 
could be significantly strengthened, and even the relevant 
structural phase transitions could occur and be observed 
under feasible experimental conditions in those metals 
whose Fermi surfaces insert nearly cylindrical belts. 

At first, we derive expressions for electron contribu- 
tions to the elastic constants. To simplify further calcu- 
lations we assume that the FS is axially symmetric, and 
the external magnetic field B is directed along the sym- 
metry axis ("z" axis of the chosen coordinate system). 
We analyze the elastic response of a metal to an exter- 
nal deformation described with the lattice displacement 
vector u(r). The effect of conduction electrons on the 
crystalline lattice arises due to a self-consistent electric 
field which appears under deformation. Also, the lattice 
deformation gives rise to an additional inhomogeneous 
magnetic field b(r), and to deformation induced correc- 
tions to the crystalline fields. Here, we omit the latter 
at first steps of our analysis to include them later. The 



emergence of the electric and magnetic fields accompa- 
nying the lattice deformation, leads to a redistribution of 
the electron density N. The local change in the electronic 
density 8N(r) equals: 



(r) = ~ac e</,(r) + aB 



N* ( ( ev>(r) + ^b(r) 



where e is the absolute value of the electron charge. 
The magnetic field b(r) satisfies the equation: 



(1) 



curlb(r) = 47rcurl<5M(r) 
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divb(r) = 0. (2) 



Here M is the magnetization vector; ( is the chemi- 
cal potential of charge carriers; ip(r) is the potential of 
the electrical field, arising due to the deformation. The 
quantity N£ included into Eq.(l) is closely related to the 
electron density of states (DOS) on the Fermi surface N^. 
The diffference between the two originates from the cor- 
relations in the electron system. Within the framework 
of the phenomcnological Fermi-liquid theory the renor- 
malizcd DOS N* has the form (See [1]): 



-V = "V^ 4^<„'(-q)r^v(q) 



E v — E v > 



(3) 
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where f v is the Fermi distribution function for quasipar- 
ticles with energies E v , and n„/„(q) is the Fourier trans- 
form of the operator of electron density in space variables. 
The renormalized operator of electron density n* vv , (— q) 
is related to the "bare" operator n vv i{— q) as follows: 



<„'(-q) = rw(-q) + 



E u 
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(4) 
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where F^f 2 are the matrix elements of the Fermi-liquid 
kernel: 

Kl> 2 = C" 2 ^^^+€" 2 (ws^ 2 ). (5) 

Here, a is the set of orbital quantum numbers, a is the 
spin number, and s is the operator of the electron spin. 

The relations (1), (2) have to be complemented by the 
condition of electrical neutrality of the system : 

SN(r) + eiVdivu(r) = 0. (6) 

It follows from the equations (1),(2),(6) that: 



curl{(l — 4-7rx)b(r)} = — 47rcurl ■ 



«9M N 



0( N 



7 divu(r) . (7) 
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FIG. 1: Schematic plot of the magnetic field dependence 
of 4:irx\\ (dashed line) and cn/c u (solid line) near a peak 
of quantum oscillations at B = Bo (T = 0). The range of 
magnetic fields corresponding to structural (Bi < B < B2) 
and/or magnetic (B[ < B < B' 2 ) instability is section- lined. 



The set of simultaneous equations (1),(2),(6) was first 
presented in previous works [1,6]. We use these equations 
to exclude b(r) and to express the potential ip(r) in terms 
of the lattice displacement vector. As a result we arrive 
at the expression for the electron force F(r) acting upon 
the lattice under its displacement by the vector u(r): 

F(r) = A b (b V(Vu(r))) + A[b x [V(Vu(r)) x b ]]. 

(8) 

Here, b is the unit vector directed along B. This result 
(8) proves that the constants Ao, Ai represent electron 
contributions to the elastic constants corresponding to 
the deformation of the lattice along the external magnetic 
field (Ao) and across this field (A). Within the adopted 
geometry these constants equal the electron terms in the 
compression elastic moduli C33 and en = c 2 2 (in Voight 
notation). On the basis of the equations (1), (2), (6) the 
expressions for these constants are derived [1,6]: 



iV 2 



Ao 

A = A [ 1 
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(9) 
(10) 



Here X|| 



dM z dM z d( 
~dB + ~diT~dB 



is the longitudinal part of 



the magnetic susceptibility; xc = 



dM z d( 



d( dB 

As follows from Eq.(9), the quantity Ao coincides with 
the compression modulus of the electron liquid. The 
structure of the quantity A is more complicated. Be- 
sides the electron compression contribution, A also con- 
tains a contribution of a different origin. This extra term 
appears due to the inhomogeneous magnetic field b(r) 
caused with the lattice deformation. This field arises due 
to the change in the magnetization of electrons caused 
by the deformation. So, the appearance of the second 
term in the expression (10) is a manifestation of a mag- 
netostriction effect. 



When the differential magnetic susceptibility enhances 
in the proximity of diamagnetic phase transition, the de- 
nominator of the second term in (10) can take on values 
close to zero. Owing to this, the quantity A significantly 
grows in magnitude. For negative A this enhancement 
brings a noticeable decrease in the elastic constant Cn 
(the latter is the sum of the "bare" elastic modulus 
and the electron contribution A). In other words, the ex- 
pression (10) reveals a possibility for softening of a lon- 
gitudinal acoustic mode propagating perpendicularly to 
the magnetic field B. This possible softening arises due 
to magnetostriction, and it results from magnetic insta- 
bility. 

When a strong magnetic field is applied to a metal, this 
gives rise to quantum oscillations in the electron DOS. 
The latter causes quantum oscillations in observables in- 
cluding magnetic susceptibility X||- At low temperatures 
(9 « 1) the magnitude of quantum oscillations increases 
so much that the oscillating term could predominate at 
peaks of oscillations. It was shown before (within the 
simple model of isotropic electron liquid) that under such 
conditions both N£ and 1 — 47TX| could go to zero near 
the oscillations peaks producing magnetic instability of 
the metal and softening of some acoustic modes [1,5,6]. 
This is illustrated in the Fig. 1. Here, magnetic fields B[ 
and B' 2 label thresholds of the magnetic instability re- 
gion, and the differential magnetic susceptibility diverges 
at these points. Singularity in the longitudinal suscep- 
tibility X|| discussed in the earlier papers [1], appears 
significantly closer to the field Bo indicating the position 
of the oscillation peak. At the same time, the structural 
instability thresholds B\ and B2 are located farther from 
Bq than B[ and B' 2 , respectively, as shown in the Figure 
1. However, these effects could be revealed in experi- 
ments only at extremely low temperatures (of the order 
of 10 mK or lower). The stringent requirements for tem- 
peratures explain why the softening of the phonon modes 
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at peaks of quantum oscillations was not observed so far. 

It is known that under certain conditions local geo- 
metrical features of the FS could significantly affect the 
electronic response of the metal. This happens when a 
dominating contribution to the response functions results 
from small "effective" segments of the FS, where "effi- 
cient" electrons are concentrated. When the effective seg- 
ments of the FS are nearly cylindrical or include locally 
flattened pieces, this gives a significant enhancement in 
the number of effective electrons, and can produce no- 
ticeable changes in the electronic response of the metal. 
The influence of locally flattened and neary cylindrical 
parts of the FS on the ultrasound attenuation rate, as 
well as on the surface impedance of a metall has been 
analyzed before (see e.g. [8-10]). Quantum oscillations 
of the electron DOS in strong magnetic fields are speci- 
fied with contributions from effective cross-sections of the 
FS. Those are cross-sections with minimum and maxi- 
mum sectional areas. Therefore the local geometry of the 
FS in the vicinities of these cross-sections has to affect 
both magnitude and shape of the oscillations [10]. This 
can influence anomalies of the clastic moduli under the 
present study and create much more favorable conditions 
for their observations in metals, as we show below. 

We consider a metal whose FS is axially symmetric in 
the vicinity of an extremal cross-section at p z = p with 
the area A ex . We assume the magnetic field B to go 
along the symmetry axis. The general expression for the 
FS curvature near p z = p could be written in the form: 



the electron DOS oscillations: 



1 d 2 A/dx 2 
2 



(11) 



Here, x — (p z — Po)/Pm] Pm is the maximum value of 
the quasimomentum component p z at the FS. Now, we 
adopt the following approximation for the cross-sectional 
area around the extremal cross-section: 



A{x)=A ex {l±b z x M ) 



(12) 



where b 2 is the dimensionless constant, (b 2 « 1) and 
the parameter I takes on values greater than 1. 

The very essence of the employed model (12) is that 
it describes an axially symmetric FS whose curvature 
turns zero at p z = pq. The expression for the curvature 
(11) showes that the approximation (12) is necessarily 
applicable to each nearly cylindrical strip on any such 
FS. Otherwise the strip has a nonzero curvature. So we 
see that the model (12) gives the general expression for 
cross-sectional area of any nearly cylindrical segment of 
an axially symmetric FS. 

Assuminig that the cyclotron quantum Ml is small 
compared to ( (j^ 1 = {Kl/Q 1 / 2 « 1) and using the 
model (12) we arrive at the following expression for the 
contribution from the nearly-cylindrical cross-section to 
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Here, tp r (8) 



sinh t-0' V 2l(b^fl l 
gamma function, and UClo is the spin splitting energy. 
Our formula (13) agrees with the results obtained for 
a precisely cylindrical FS/and a Fermi circle in a two- 
dimensional conductor (see e.g. [11,12]). We arrive at 
the corresponding results within the limit I — > oo. In 
general case we can treat "Z" as a phenomenological pa- 
rameter included in the model (12). Actual values of "Z" 
could be discovered in experiments where the FS local 
geometry is revealed. This is the only trustworthy way 
to estimate the above parameter for a particular metal. 
First principle calculations are not accurate enough to 
produce reliable results on fine geometrical features of 
FSs although their shapes as a whole are well known. 

Within the isotropic model the cross sectional area 
is described with the expression (12) where I = b 2 = 
1, po = 0, and p m is the radius of the Fermi sphere in 
quasimomemta space. In this particular case the oscil- 
lating function A takes on a well-known form: 



r(i/2Z) 
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(14) 

Oscillations described by (13) and (14) differ in phase 
as well as amplitude. The amplitude of usual oscillations 
given by (14) is of the order of 7~ 1 0~ 1 / 2 , while (13) gives 
a magnitude of the order of -y -1 /^ 1 - 2 ')/ 2 ' . Therefore, 
the amplitude of oscillations related to the extremal sec- 
tion of zero curvature is approximately (^j~ 1 9 1 / 2 )(- 1 ^ l '>/ 1 
times greater than that of the usual quantum oscillations. 
As a result, the contribution due to an extremal section 
with zero curvature can be considerably (more that ten- 
fold) greater than contributions due to other extremal 
sections, and the function A can reach values of the or- 
der of unity at peaks of oscillations even at 9 <~ 1. On 
these grounds we conclude that the most favorable con- 
ditions for observation of softening of elastic moduli at 
peaks of quantum oscillations occur in metals whose FSs 
include nearly cylindrical segments. We consider such 
FSs in further analysis. 

Assuming 7 << 1 we replace matrix elements included 
in the Fermi-liquid kernel with their semiclassical analogs 
<p(p, p') and V(P)P') which depend on quasimomenta of 
interacting conduction electrons p and p'. For axially 
symmetric FSs the Fermi-liquid functions take the form: 

<P(P,P') = Voo+PzP' z Voi + (P±p'±)( i Pw+PzPWu) : > (I 5 ) 



^(P,P') = i>oo+PzP' z ipoi + (p±P±)(i>w+PzP' z i>ii)- (16) 
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To make our analysis more thorough we include defor- 
mation terms in the force (8) exerted by the conduction 
electrons on the lattice. We present the components of 
the deformation potential, as: 



A Q /?(p) = A Q/3 + GU a/3 (p), 



(17) 



where A a p is a tensor whose elements do not depend on 
p, G is a dimensionless constant, and Tl a p is the tensor 
of the electron momentum flux density. 

The adopted approximations lead to the following ex- 
pressions for the corrections to the elastic constants: 

- - ^ 2 [1 + G/(1 + A 2 )] 2 A 

C11 = C22 = ~T l + (l + ^-4W)A' ( 8) 



Cl2 



N 



G 



g \1 + A 2 J 1 + (1 + T^-4ttxo7 4 )A 



• (19) 



where g is the electron DOS in the absence of the external 
magnetic field, 



C33 = - 



N 2 [1 + G(1 + ^ 2 ) 2 ]A 
g 1 + (1 + 1F)A 



(20) 



Xo is related to the Landau diamagnetic susceptibility 
(the latter equals — ^Xo); an d the constant W originates 
from the Fermi-liquid interaction: 



1 + B Q 1 + B 2 1 + A 2 



(21) 



Here, dimensionless coefficients B ,B 2 and A 2 are re- 
lated to the Fermi-liquid parameters V'oo, V'lii an d <fin as 
follows: 



B ° = (2W / MPz)dPz 



B- 



ipuS(p z )p z dp z 



(2tt/j) 3 

A 2 = j^j-^ J ViiS(p z )p z d Pz . (22) 

It is noteworthy to mention that the shear modulus c\ 2 
is also affected due to magnetostriction. This was not 
discovered in the early analysis of [1] where the deforma- 
tion contributions into the electron force were omitted. 
As for elastic constant C33, it also can be reduced at the 
peaks of low temperature quantum oscillations but this 
effect is not related to magnetic instability. The expres- 
sion for £33 (20) does not include the contribution arising 
due to magnetostriction. The possible softening of C33 is 
directly connected with the behavior of the electron DOS 
under strong magnetic fields at low temperatures. Such 
instability was predicted before [7]. We have to remark 
here that the interaction among electrons significantly 
influences all above effects. The value of the constant 
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FIG. 2: Magnetic feld dependencies of the elastic constant 
Cn near the diamagnetic phase transitions. The curves plot- 
ted for: = 1, the parameter I takes on values 8,4,2 from 
the left to the right (left panel); I = 4, 6 takes on values 
1,2,3 from the left to the right (right panel). For all curves 
N = 10 21 cm~ 3 , 7 2 = 10 3 , Bo = 10T. 



W which accumulates effects of electron-electron inter- 
actions within the framework of the Fermi liquid theory 
could significantly influence the temperature range where 
both magnetic and lattice instabilities occur. 

The function A describing quantum oscillations is 
given by Eq. (13) and the amplitude of oscillations 
may become comparable with unity at moderately low 
temperatures provided that the FS shape reveals a fair 
proximity to a cylinder near the extremal cross sec- 
tion. For example, at I = 3, hfl/( <~ 10~ 3 , and 
B ~ 10T, the condition 7 -V*0(i-2O/2l _ 1 could be 
satisfied at temperatures of the order of 1 K. In usual 
metals N <~ 10 21 -j- 10 22 cm -3 , and the term 47rxo7 4 m 
the denominator of (18) can take on values of the or- 
der of 10 at 7 2 <~ 10 3 . Also, we introduce the quantity 
L 2 = N 2 (l + G/(l + A 2 )) 2 /gc°, where c° is the relevant 
"bare" elastic constant. The ratio N 2 /gc° in typical met- 
als is rather small (N 2 /gc° ~ 10~ 2 10 _1 ). However, 
the values taken on by L 2 could be noticeably greater 
than that due to the deformation constant G. The latter 
mostly accepts values of the order of unity. So, we have 
grounds to expect the values of L 2 to be of the order of 
10 _1 rather than 10~ 2 . 

To proceed in the analysis of the experimental feasi- 
bility of the effect, we numerically evaluate the decrease 
in the elastic constant Cn using our result (18) and the 
above estimations of the parameters included there. The 
results are shown in the Fig. 2. We see that close enough 
proximity in the shape of an effective strip on the FS to a 
cylinder gives rise to the structural instability near peaks 
of oscillations at 8 = 1. Also, it is demonstrated that the 
effect is washed out due to the further rise in tempera- 
ture. 

Electronic contributions to the velocity of ultrasound 
waves propagating in metals are simply related to the 
elastic constants. It follows from our results that the 
velocity of a longitudinal sound could depend on the di- 
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rection of its propagation. In the vicinity of the Condon 
instability the velocity of sound propagating perpendicu- 
larly to the magnetic field B could be noticeably reduced 
compared to the velocity of sound propagating along B. 
Again, we may expect this effect to appear in metals 
whose FSs include nearly cylindrical segments. 

Finally, it is known that the effect of magnetostric- 
tion can cause softening of some phonon modes in met- 
als near the Condon magnetic instability. This effect can 
appear even in an isotropic metal [1,6]. The point of 
the present work is that the effects could be significantly 
strengthened when the immediate vicinities of some ex- 
tremal cross-sections of the FS are nearly cylindrical in 
shape, so that the FS curvature turns zero at these cross- 
sections [13]. there is an Real metals mostly have non- 
spherical and complicated in shape FSs. At present the 
main geometric characteristics of the FSs, such as their 
connectivity, are there is an well studied. On the con- 
trary, local geometric features of the FSs has not been 
investigated in detail so far. However, there is an exper- 
imental evidence that "necks" connecting quasispherical 
pieces of the FS of copper include nearly cylindrical belts 
[11]. When the magnetic field is directed along the axis 
of a "neck" (for instance, along the [111] direction in the 
quasimomenta space), the extremal cross section of the 
"neck" could be expected to run along the nearly cylin- 
drical strip where the FS curvature turns zero. It is also 
likely that the FS of gold possesses the same geometrical 
features for it closely resembles that of copper. Another 
kind of materials where we can expect the low tempera- 
ture softening of phonon modes to be manifested includes 
some layered structures with metallic-type conductivity 
(e.g. a-(BEDT-TTF) 2 MHg(SCN) 4 group of organic 
metals). Fermi surfaces of these materials are sets of rip- 
pled cylinders, isolated or connected by links [14]. Based 
on the experiments on cyclotron resonance in these or- 
ganic metals [15], it was shown that the cylinders could 
have nearly cylindrical strips [16]. For all above listed 
substances we can expect the effect to be revealed at rea- 
sonably low temperatures (<~ IK) and reasonably strong 
magnetic fields (1 -j- 10T). 

The effect is expected to be very sensitive to the ge- 
ometry of an experiment for extremal cross sections of 
the FS run along nearly cylindrical belts (if any) only 
at certain directions of the magnetic field. When the 
magnetic field is tilted away from such direction, the 
extremal cross section slips from the nearly cylindrical 
piece of the FS. This does not cancel possibilities of the 
effect in principle but makes requirements for tempera- 
tures dramatically stringent, as it was discussed above. 
The present analysis is carried out within the model of 
axially symmetric FS. The obtained results could be ap- 
plied to actual metals when the magnetic field is directed 
in parallel with a high order axis of symmetry of the crys- 
talline lattice of a metal. Otherwise, it is very difficult to 
separate out electron contributions to particular elastic 



constants. In summary, the effect of softening of phonon 
modes at peaks of quantum oscillations in metals could 
be expected to be observed in experiments. The most 
favorable conditions for the effect exist in metals whose 
FSs include nearly cylindrical pieces. The efffect could 
be revealed in these metals for some particular directions 
of the external magnetic field providing that an extremal 
cross section belongs to a quasicylindrical strip on the 
FS. 
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